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ABSTRACT 

Context. The Cygnus region harbours a huge complex of massive stars at a distance of 1.0-2.0kpc from us. About 170 O stars are 
distributed over several OB associations, among which the Cyg 0B2 cluster is by far the most important with about 100-120 O stars. 
These massive stars inject large quantities of radioactive nuclei into the interstellar medium, such as ^*A1 and ^''Fe, and their gamma- 
ray line decay signals can provide insight into the physics of massive stars and core-collapse supernovae. 

Aims. Past studies of the nucleosynthesis activity of Cygnus have concluded that the level of ^^Al decay emission as deduced from 
CGRO/COMPTEL observations was a factor 2-3 above the predictions based on the theoretical yields available at that time and on 
the observed stellar content of the Cygnus region. We reevaluate the situation from new measurements of the gamma-ray decay fluxes 
with INTEGRAL/SPl (presented in a previous paper) and new predictions based on recently improved stellar models. 
Methods. We built a grid of nucleosynthesis yields from recent models of massive stars. Compared to previous works, our data in- 
clude some of the effects of stellar rotation for the higher mass stars and a coherent estimate of the contribution from SNlb/c. We 
then developed a population synthesis code to predict the nucleosynthesis activity and corresponding decay fluxes of a given stellar 
population of massive stars. 

Results. The observed decay fluxes from the Cygnus complex are found to be consistent with the values predicted by population 
synthesis at solar metallicity; and yet, when extrapolated to the possible subsolar metallicity of the Cygnus complex, our predictions 
fail to account for the INTEGRAL/SPl measurements. The observed extent of the 1809 keV emission from Cygnus is found to be 
consistent with the result of a numerical simulation of the diffusion of Al inside the superbubble blown by Cyg 0B2. 
Conclusions. Our work indicates that the past dilemma regarding the gamma-ray line emission from Cygnus resulted from an over- 
estimate of the 1809keV flux of the Cygnus complex, combined with an underestimate of the nucleosynthesis yields. Our results 
illustrate the importance of stellar rotation and SNlb/c in the nucleosynthesis of ^*A1 and ^"Fe. The effects of binarity and metallicity 
may also be necessary to account for the observations satisfactorily. 

Key words. Gamma rays: observations - open clusters and associations: individual: Cyg OB - Stars:early-type - ISM:bubbles - 
Nuclear reactions, nucleosynthesis, abundances 



1. Introduction 

Conventionally defined as stars with initial masses above 8-10 
, and up to about 120 Mq, massive stars play a fundamental role 
• in astrophysics because of the very powerful behaviour they ex- 
' hibit during their lives, from the early main-sequence to their 
final explosion (and even beyond thanks to the compact objects 
they leave behind). Usually clustered in OB associations of a 
few tens to several thousand members, massive stars severely 
. affect their environment on the hundred-parsec scale through 
their sustained mechanical and radiative energy output: strong 
luminosity in the Lyman continuum, very energetic winds and 
mass loss episodes, and supernova explosions. In addition to me- 
chanical energy, these stellar outflows also carry the products of 
nuclear burning and thus enrich the interstellar medium (ISM) 
with heavy elements. The physical and chemical evolution of 
the Galaxy is driven to a considerable extent by massive stars 
and a wide variety of Galactic objects follow from their evolu- 
tion: supernovae (SNe) and their remnants (SNRs), neutron stars 
or stellar-mass black holes, some of which manifest in the form 
of pulsars or magnetars. X-ray binaries, or microquasars. 

A thorough understanding of the processes that rule the evo- 
lution of massive stars and their death as core-collapse SNe is 



therefore inevitable in astrophysics. Amongst the various obser- 
vational approaches employed to explore the physics of massive 
stars and their SNe, gamma-ray line astronomy can be a very 
valuable tool. By studying the characteristic decay radiation of 
certain radio-nuclides produced in the hydrostatic and/or explo- 
sive burning stages of massive stars, gamma-ray line astronomy 
provides indirect access to the physical conditions that governed 
their synthesis and ejection in the ISM. Among the few radio- 
isotopes accessible to the present-day gamma-ray telescopes, the 
long-lived Al and ''"Fe hold potential for interesting constraints 
on the evolution of massive stars. 

The Galactic decay emission from Al was extensively stud- 
ied in the 1990s with the CGRO/COMPTEL instrument and the 
resultant cartography of the ^''Al decay signal at 1809keV re- 
vealed that most of the emission comes fro m the inner Galaxy 
and from a few ne a rby star-forming regi ons jDiehl et alj|1994t 
Idel Rio et all 1 19961: IPmichke et aLlfeoci ll) . The intensity distri- 
bution observed with COMFTEL allowed it to be established 
that massive stars are the domina nt source of ^^AI in the Galaxy 
(Prantzos & Diehl 1996; Knodl seder et all 1 19991) . In the early 
2000s, the first detection of the Galactic ^''Fe decay signal at 
1173 and 1332 keV was achieved from RHESSI observations 
(■Smith.2004.) . with a measured flux of about 15% of the ^""Al de- 
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cay flux. Although the decay emission was too weak to achieve 
a detailed cartography and thus confirm that massive stars are 
the main source of ^"Fe, the combined analysis of Al and ''"Fe 
decay signals appeared as a potentially strong constraint on the 
processes that govern the nucleosynthesis in massive stars and 
SNe. 

The launch in 2002 of the INTEGRAL gamma-ray obser- 
vatory, equipped with the SPI high-resolution spectrometer, has 
opened the possibility of extending the analysis of the gamma- 
ray line emission from the Galaxy and add spectrometric in- 
formation to the COMPTEL results. Both the ^^Al and *"Fe 
decay signals have now been unambiguously detected by SPI 
on the Galactic scale (Diehl et al. 2006; Wang et al.,,2007 ). The 
stronger ^^Al emission made it possible to analyse its morphol- 
ogy in more detail and appeared particularly prominent in some 
specific star-forming regions of the Galaxy, in agr eement with 
the COMPTEL observations. In a previous work ( Martin et al.l 
12009, hereafter Paper I), we showed that the Cygnus region ex- 
hibits a relatively strong signal of Al decay, so we investigated 
further the gamma-ray line emission of both Al and ^"Fe from 
Cygnus. 

The Cygnus region and its peculiarities were presented in 
detail in iPaper J and we only repeat the main characteristics of 
the region here and the principal arguments for studying it in 
nuclear gamma-ray lines. The Cygnus region, which we con- 
sider to be that part of the Galactic disk located between lon- 
gitudes 70° and 90°, harbours a very impressive concentration 
of massive stars at o nly ~ I-2kpc from us. A bout 170 O stars 
have been recorded (iKnodlseder et al.l |2002|) . distributed in 6 
OB associations and a dozen open clusters, among which the 
spectacular Cyg OB2 cluster tha t alone com prises 100-120 O 
stars dKnodlseder 2000; Comero n et al. I I2OOI . In the following, 
we collectively refer to these OB associations and open clusters 
(the details of which are presented later) as the "Cygnus com- 
plex". With this term, we aim at gathering all these rich and 
nearby stellar groups that are relatively correlated in terms of 
distance (between 1 and 2 kpc) and thus contribute to the speci- 
ficity of the Cygnus region. This "Cygnus complex" appears to 
be quite young, as shown by the lack of radio supernova rem- 
nants (SNRs) or pulsars. The Green catalogu43 gives about 10 
radio SNRs betwee n 70° and 90° , but a ccording to a compila- 
tion of distances by iKaplan et al ] (12004 . none of these objects 
can be related to the Cygnus complex. IWendker et alJ (Il99lh 
concluded from the sensitivity of their radio survey of the re- 
gion that no other radio SNR exists up to a distance of 10 kpc. 
Additional evidence for the relative youth of the complex comes 
from the absence of radio pulsars tow ards Cygnus up to 4 kpc , 
as indicated by the ATNF catalogueQ ([Manchester et al.ll20O5l) . 
although recent observations with the LAT telescope aboard the 
Fermi Gamma-Ray Space Telescope have led t o the discovery of 
the young energetic pulsar PSR J2032-I-4127 (Abdo et al.ll2009l) 
that is likely associated with Cyg OB2 (Camilo et al. 2009). For 
those reasons, the Cygnus region appears well-suited to studying 
massive stars in their hydrostatic phases. 

The Cygnus population of massive stars is close and rich 
enough to ensure sufficient gamma-ray line fluxes and there- 
fore constitutes a good opportunity to test our understand- 
ing of the structure and evolution of massive stars, as re- 
vealed by their nucleosynthetic activity. A number of stud- 



' Green D.A., 2006, "A Catalogue of Galactic Supernova Remnants 
(2006 April version)". Astrophysics Group, Cavendish Laboratory, 
Cambridge, UK ( http://www.mrao.cam.ac. uk/s urveys/snrs/| . 
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ies have been devoted to the nucleosynthesis of ^^Al and 
''"Fe in the Cygnus region (see for example del Rio et aL 
1996 ; Pluschke et al. 2000; Cervino et al. 2000; Pluschke et at 
12002! : iKnodlseder et al J 12002 ) and the latest studies agree that 
the level of the 1809keV emission from Cygnus as deduced 
from CGRO/COMPTEL observations is a factor of 2-3 above 
the predictions, when the latter are based o n the theore tical 
viel ds bv Meynet et al . ( 1997), Woosley & W eaveii ( Il995h and 
Wo oslev et al.l d 1 9951) and on the observed stellar content of the 
Cygnus region (includ ing the upward rev i sion of the stellar pop- 
ulation of Cyg OB2 bv lKnodlsedeil2000h . lPluschke et al.1 (l2000h 
show that the discrepancy could be alleviated if one includes 
the enhanced yields from massive close binary systems in the 
calculation. Alternatively, iPliischke et al.l ( |2002|) argue that the 
problem may well come from an underestimate of the stellar 
content of the Cygnus OB associations due to a strong inter- 
stellar e xtinction, as was shown to be the case for Cyg OB2. In 
iKnodlseder et alj (12002 ). the authors speculate that future stellar 
models, especially those that include the efl'ects of stellar rota- 
tion, could solve the dilemma. 

In Paper I, we used the INTEGRAL/SPI observations to re- 
assess the ^''Al and ^'^Fe decay emission from the Cygnus re- 
gion. In particular, we disentangled in a coherent way the con- 
tribution of the nearby clustered OB population of the Cygnus 
complex from the contribution of the more diffuse and spa- 
tially extended non-clustered population, which is likely to con- 
tribute as well to the observed emission. We now compare our 
INTEGRAL/SPI observations of the gamma-ray line emission 
of the Cygnus complex with theoretical predictions based on the 
latest models of stellar nucleosynthesis. In recent years, stellar 
models have been improved with modifications, such as the in- 
clusion of stellar rotation or the revision of the mass loss rates 
and nuclear cross-sections. These upgrades have already proven 
to bring simulations closer to observations. For instance, mod- 
els including rotation can account for changes in surface abun- 
dances oc curring during t he main sequence phases ( H unter et al] 
2008; M aeder et al.ll2009l) and for the observed variations of the 
number of Wolf-Rayet to O-type stars as a function of metallic- 
ity in constant star formation regions (Mevnet & Maeder 2005]). 
In the following, we show that the new stellar models, combined 
with the work done on observations in Pa per L considerably al- 
leviate the long standing discrepancy between theoretical and 
observed decay fluxes from Cygnus. 

2. Observational constraints from INTEGRAL/SPI 

In iPaper J, we presented the analysis of about 4 years of 
INTEGRAL/SPI observations amounting to a total effective ex- 
posure time of 63.2 Ms. Most of this exposure is confined to the 
Galactic plane, with 10.8 Ms specifically covering the Cygnus 
region. In this section, we review the main results that were ob- 
tained on the ^*A1 and ''"Fe decay emission. 

2.1. The 1 809 keV emission from ^^Al 

The ^^Al isotope is produced predominantly by massive stars 
at various stages of their existence and it is released in the 
ISM by both the stellar winds (of Wolf-Rayet stars espe- 
cially) and the SNe. Nucleosynthesis of ^^Al follows from the 
^^Mg(p,7)^^Al reaction, which mostly occurs during H central- 
burning (with ^^Mg coming from initial metallicity), C/Ne shell- 
burning (with ^^Mg coming from the CNO initial metallic- 
ity through the ^'^n{a,yy^¥{&-Y'^0{a,yf'^Ne{a,nf^Mg chain) 
and explosive burning of the C-shell (with ^^Mg coming from 
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^'*Mg(n,y)^^Mg, ^"^Mg being a product of C/Ne burning). 

The mean ^^Al lifetime of about 1 Myr allows it to diffuse 
in the local ISM before decaying into ^^Mg. As a consequence, 
the corresponding 1809keV emission is of diffuse nature with 
an overall distribution that follows that of massive stars in our 
Galaxy, as also traced for instance by the microwave free-free 
emission from HII regions or the infrared emission from heated 
dust (Knodlsederet al. 1999). Most of the Galactic 1809 keV 
flux is thus confined to the Galactic plane and inside the inner 
regions, except f or a nota ble emission from the Cygnus and Sco- 
Cen regions (see lPaperl) . 

From the INTEGRAL/SPI data, we have found that the dif- 
fuse 1809keV emission from the Cygnus region is well repre- 
sented by a 3° X 3° 2D Gaussian intensity distribution (which 
implies a characteristic angular size of 9-10°) centred approxi- 
mately on the position of the Cyg OB2 cluster This confirms ex- 
pectations that the latter dominates the energetics and nucleosyn- 
thetic activity in this area. It should be emphasised, however, 
that we only have weak constraints on the maximal extent of the 
emission due to intrinsic limitations of the coded mask imaging 
system of the SPI instrument and yet, the size we found is consis- 
tent with the results obtained from the CGRO/COMPTEL obser- 
vations (iPliischke et al...2001) . The total flux from the extended 
source in the Cygnus direction is (6.0 + 1 .0) x 10"^ ph cm"^ s ' , 
where (3.9+ 1.1) x lO^^ph cm ^ s ' is attributed to the Cygnus 
complex and the remaining flux stems from Galactic disk fore- 
ground and background emissions. 

Thanks to the high spectral resolution of the SPI instrument, 
we were able to study the profile of the 1809keV line. We ob- 
served a line position corresponding to a radial velocity of 33 
+ 66 km s"' , which agrees with the overall radial motions in the 
direction and at the distance of the Cygnus complex, as indi- 
cated for ins tance by spectrometric observations of the CO line 
dPame et alJ feOOl). We also measured a slight line broadening 
that, within the statistical uncertainty, is consistent with ^*A1 be- 
ing still or in a medium with turbulent or expansion velocities of 
up to 100-200 km s-i. 

2.2. The 1 1 73 and 1332keV emission from ^°Fe 

The massive star origin of the ^''Fe isotope is far less established 
than for ^*A1. Th e nucleosynthesis of *"Fe by massive stars 
was addressed by Timmes et al. ( 1995) and extensively reviewed 
by (Limongi & Chieffi 2006). From these works, we know that 
massive stars are a source of ^"Fe through the ^^Fe(n,7')^'^Fe re- 
action that occurs during He shell-burning, C shell-burning, and 
explosive burning of the C-shell. In these processes, the seed 
amounts of Fe come from the initial metallicity and so do the 
neutrons, through the ^^Ne(Q',n)^''Mg reaction (with ^^Ne com- 
ing from the initial CNO, see 12.11 ). The ^"Fe thus produced by 
massive stars is then released in the ISM only through SN ex- 
plosions and no wi nd contribution is expect ed. Other objects 
suc h as AGB stars (Karakas & Lattanzio 2007) or SNe of type 
la dlwamoto et al.l il999) are also potential producers of ^*^Fe. 
Although the relative contribution of each class of objects re- 
mains entirely unconstrained on the observational level, the 
available models for the various potential sources indicate that 
core-collapse SNe should provide most of the Galactic ^"Fe. 

The observed intensity of the Galactic decay emission at 
1173 and 1332keV i s a factor of 7-8 lower flian the 1809 keV 
emission from ^^Al ( Harris et al.1 120051: IWang et al.ir200 7'. and 
[Paper I), so the detection of any ^''Fe signal from Cygnus may 
at first sight be difficult given the sensitivity of INTEGRAL/SPI. 
In addition, the Cygnus complex seems quite young, and it may 



well be that its stellar population has so far not produced a single 
SN, which would mean no release of ^"Fe into the ISM, hence 
no 1173/1332keV signal from the Cygnus complex at all. As a 
matter of fact, we did not observe any 1173/1332keV emission 
from the Cygnus direction and derived a 2cr upper limit on the 
''"Fe decay flux o f 1.6 x lO'^phcm '^ s"' (consistent with the 
value obtained bv lWang et al.ll200'7l) . 

3. Population synthesis 

In this section we present the most important aspects of the 
population synthesis simulations we performed to compare our 
INTEGRAL/SPI observations with. In this eff'ort, we took ad- 
vantage of the recent progress by two research groups in the field 
of massive stars nucleosynthesis. In the following, we introduce 
the specificities of each grid of stellar models and provide some 
detail about their implementation in our predictive tool. 

3. 1 . The stellar models 

Once the massive stars had been recognised as the dom- 
inant source of the Galactic ^^Al (Prantzos & Diehl 199^ 
iKnodlseder et aljl999h . the question of its origin shifted to iden- 
tifying the subset of massive stars and/or evolutionary phases re- 
sponsible for its ejection. For historical reasons, mostly linked 
to the evolution of the stellar models, two "competing" sources 
were discussed: stars less massive than about 35 M© through 
their SNII explosions, and stars more m assive than 35 Mp 
through their WR winds (see the review by iLimongi & Chieffil 
2006, hereafter LC06). The debate then settled on the relative 
contribution of SNII and WR to the Galactic ^^Al budget and, in 
this situation, the ^"Fe appeared as a promising discriminating 
agent since it is released only by SN explosions. 

The work oflEcoe' has broadened the issue thanks to a new 
grid of stellar models. The authors computed the nucleosynthe- 
sis yields for a wide range of non-rotating solar metallicity mod- 
els, with initial stellar masses from 11 to 120 M©. The models 
were followed over their entire life, including both the hydro- 
static phases and the explosive burning episodes of the super- 
nova explosion. The main outcomes of the LC06 study are: 

1 . Over the mass interval considered, most of the ejected ^^Al 
has explosive origins. For the particular case of the most 
massive stars, typically the WR progenitors, the production 
of ^^Al due to explosive burning of the C/Ne shells exceeds 
the quantity synthesised by central H-burning and carried 
away by the stellar winds. 

2. The production of ^"Fe (by C/Ne convective shells burning 
mostly) of the most massive stars that explode as SNIb/c is 
up to a factor 10 higher than the yields associated with the 
lower mass stars that explode as SNII. 

It has always been clear that a dichotomy between WR and SNII 
as the sources of ^*A1 is inadequate and that the whole contribu- 
tion of each massive star should be accounted for to understand 
the origin of both the ^^Al and ^"Fe. Moreover, the initial metal- 
licity may have a strong impact on the relative contributions of 
WR winds and SN explosions to the production of ^^Al, as will 
be illustrated later. 

In parallel to the work of lLC06l grids of stellar models 
including stellar rotation have revealed that the latter strongly 
affect s the evolut i on and nuc leosynthesis of massive stars. 
Me vnet & Maedeil (l2003l l2005h have computed the evolution of 
rotating models of stars with initial masses from 20 to 120Mo 
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Fig. 1. Supernova yields of A l and Fe as a function of the CO 
core mass (from the models of lLC06l) . The se curves are used to 
associate the purely hydrostatic models of iMMPOSi with super- 
nova yields. 



up to the end of central He-burning. These rotating models were 
computed for velocities compatible with the mean observed ve- 
locities for OB stars, and the impact o f stellar rotation on the pro- 
duction of ^''Al has been examined bv lPalacios et al.l (l2005h . The 
ma in outco mes of these works (hereafter collectively referred to 
as are: 

1. Rotation increases the convection and the circulation of 
chemical species inside the star. This implies a larger initial 
reservoir of ^^Mg available for proton capture (in the case 
of non-zero initial metallicity) and an earlier ejection of the 

Al produced by central H-burning. 

2. Rotation increases the mass loss by altering the surface grav- 
ity, the effective temperature, and the opacity of the envelope. 
This, combined with the more efficient diffusion of the burn- 
ing products to the stellar surface, results in an earlier transi- 
tion to the WR phase and a decrease in the minimum initial 
mass to become WR. As a result, the ^^Al contained in the 
He core (following central H burning) is released sooner and 
by more stars. 

In summary, stellar rotation alters the ^''Al yields in two ways: 
the quantities ejected in the WR winds are increased and the time 
profile of the injection into the ISM is different as well. While 
the details of the ^*A1 release with time are not important for a 
large, stationary system such as the Galaxy, for which only the 
total yields per star are relevant, the Cygnus complex of massive 
stars is a young, evolving region and the comparison of observa- 
tions with predictions therefore includes a temporal dimension. 
Therefore, a reliable population synthesis model should not only 
be able to reproduce the observed 1809keV flux from Cygnus 
but also reach an agreement within a time range that is consis- 
tent with the age estimates of the specific stellar clusters. 

Besides the above improvements, both grids of stellar mod- 
els are also based on revised mass loss rates taking the effect of 
clumping in the wind s into account and on revised nuclear re- 
action rates. MMP05 also explored the influence of the initial 
metallicity Z by comp uting fo ur grids at Z = 0.004, 0.008, 0.02, 
and 0.04. The grid of lEcOel was only computed for a Z=0.02 
solar metallicity. 



3.2. The grid of yields 

The lLCOd grid is complete and homogeneous but does not in- 
clude stellar rotation and is only available for solar metallicity, 
while the MMP05 grid includes the effects of stellar rotation 
and metallicity, but only covers the longest hydrostatic phases 
of those stars initially more massive than 20 M©. To make the 
most realistic predictions of the nucleosynthesis activity of the 
Cygnus complex, we therefore built our own grid of yields from 
the works of LC06 and MMP05. 

Before entering the description of our grid, we emphasise 
the terminology that will be used thereafter. We divided the nu- 
cleosynthesis contributions into "wind yields" and "SN yields", 
which is not equivalent to a separation into hydrostatic and ex- 
plosive yields. Wind yields pertain to ^*A1 alone and are obvi- 
ously of hydrostatic origin, whereas SN yields concern both iso- 
topes and have both hydrostatic and explosive origins. For in- 
stance, the ^**Fe released in a SN explosion comes from the hy- 
drostatic burning of the He and C convective shells but also from 
the explosive burning of the C shell. We also note that the wind 
yield is a progressive, time-dependent outflow of ^*A1, whereas 
the SN yield is considered an instantaneous release of ^*A1 and 
«'Fe. 

The grid we used in our population synthesis code is based 
on the following assumptions: 

1 . The lower mass s tars, fr om 1 1 up to 20 M©, are represented 
by the models of lLC06l for both their wind and SN contri- 
butionfl It has also been assumed that these LC06 models 
are valid at metallicities other than solar. We will discuss this 
assumption later on. 

2. The higher mass stars, from 20 to 120 Mq (typically those 
that become WR stars), are represented for their wind contri- 
bution by the rotating models of MMP05 for the four metal- 
licities considered by these authors. Each of these stellar 
tracks is then associated to one of the SN yields of lLCOel 
through its CO core mass (as explained below). 

In the composition of our grid, we assumed that the SN yield is 
strongly linked to th e final CO core mas s of the model (a hy- 
pothesis also used bv lCervino et al. I I2OOOI) . This can first be jus- 
tified by the fact that, at the moment of core-collapse, most of 
the ^*A1 and ^"Fe synthesised so far is locked in the CO core as 
a result of the He/C/Ne shell-burning episodes. Second, it seems 
reasonable to consider that the amounts of ^^Al and ^''Fe pro- 
duced explosively (through explosive burning of the C-shell) de- 
pend on the size of the CO core (even if other parameters, such 
as the chemical profile and the mass-radius relation, may be im- 
portant). From the LC06 grid of models, it is possible to draw 
relations between the SN yields and the CO core masses (see 
Fig ID- These relati ons are t hen used to associate the purely hy- 
drostatic models of MMP05 with SN yields through their final 
CO core ma ss that, for same initial stellar mass, differs from that 
obtained bv iLCOd . Indeed, for stars below about 60 M©, stellar 
rotation increases the final CO core mass because of an enhanced 
convection; above 60 Mq, rotation causes a decrease in the CO 
core mass because of an enhanced mass loss (iMevnet & Maede3 
2003. 2005). 

This approach takes into account how stellar rotation can af- 
fect the SN yields through its impact on the stellar evolution. It 
should still be emphasised that the SN yields used in our grid 



^ In practice, however, the -''Al wind contribution of these stars is 
negligible because they experience little mass loss compared to the 
higher mass stars. 
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Fig. 3. 2^A1 supernova y ields fro m |LC06| and lWHOTl compared 
to the o lder results fr om lwW95l . Also plotted are the yields for 
MM97' and 'MMP05' extrapolated through CO core mass from 
WLW95 and LC06. respectively. 



come only from non-rotating models. The impact of stellar ro- 
tation on the nucleosynthesis of the latest evolutionary stages 
remains unexplored up to now. ILC06 pointed out how important 
the shell convection efficiency is to the synthesis of ^''Al and 
^"^Fe, but how stellar rotation affects the nucleosynthesis of the 
He/C/Ne shell-burning episodes is still an open question. Our 
limited knowledge on that point therefore constitutes a potential 

source of error for our SN yields. 

In F igs. |2] and [3] we show th e ^''Al yie lds for lLC06l and 
IMMP05I (where the SN yields for iMMPOSi were extrapolated 
through CO core mass, as explained above), compared to the 
previous generation of yields from Mevnetetal. (1997, here- 
after ]VIM97), Woosle v & Weaveii (I1995L hereafter WW95) and 
IWooslev et all (i1995l hereafter WLW95). This previous gener- 
ation of yield s was u sed i n a population sy n thesis analysis by 
ICervino et aP (|2000') and ' Kn5dlseder et alJ (|2002:), who con- 
cluded that they could not account for the 1809keV flux from 
Cygri us observed by CGRO/COMPTEL (see also Pluschke et al. 
I200lh . In Fig.H we show the «'Fe yields foriLC06. and ,MMP05i 
(where again the SN yields for IMMP05I were extrapolated 
through CO core mass). 

In all t hree figures, we a l so com pare the set of yields ob- 
tained by IWooslev & Hege^ 120071 hereafter WH07). These 
yields come from a large grid of stellar models (from 13 to 
120 M0) that basically includes revised mass loss rates and re- 
vised nuclear reaction rates, but no rotation. The authors used 
an initial metallicit y of Z=0 . 016, c orresponding to the revised 
solar abundances of lLodderi (l2003l) . and computed both hydro- 
static and explosive phases. In its general features, the |WH07| 
grid is very similar to the LC06 grid, except for the different 
initial metallicity. In the details, however, many aspects of stel- 
lar physics were implemented differently (such as the mass loss 
prescriptions, see bel ow). As seen in the next paragraph s, the dif- 
ferences between the IWH07 ' yields and the LC06 and' MMP05l 
yields illustrate the extent to which the implementation of the 
stellar physics can affect the nucleosyn thesis results. 

In Fig. 121 we see that the IMMP05 yields obtained from ro- 
tating models are more enhanced than the yields of MM97, and 
this despite that the mass loss rates used in the MMP05 models 
are lower than in the MM97 models. This illustrates the impor- 
tance of the effects of rotation that, as explained in 13.11 favours 
the production of ^''Al and its subsequent ejection by the stel- 



lar w inds. The |WH07| wind yields are all above the yields of 
IlC06 (up to a fac tor of 5 in some cases) but always remain be- 
low the IMMP05I yields. This fac t very likely comes from the 
absence of rotation in the IWH07I models, but the difference be- 
tween WH07 and lLC06l could seem surprising. Indeed, a lower 
initial metallicity like Z=0.016 means that all other parameters 
being unchanged, a reduced initial reservoir of ^^Mg available 
for proton capture, hence low er ^^Al wind yields, but a thorough 
comparison ofthe WH07'and|LC0l models is beyond the scope 
of this paper. The same rates were used for the main rea ctions 
governing Al production, but the higher mass models of WH07I 
were found to experience more mass loss (as discussed in more 
detail below). This may explain the situation above 40 M©. For 
the lower-mass models, however, the higher ^''Al wind yields 
obtained with a lower initial metallicity may stem from the dif- 
ferent treatment of the convective physics and/or to the different 
mass loss rate prescription during main sequence or red super- 
giant stages. 

In Fig. |3] the ^^Al SN yields of LC06' are plotted together 
with the SN yields extrapolated through the CO core mass for 
the MMP05 rotating WR models. One clearly sees the effect of 
rotation on the size of the CO core, hence on the extrapolated SN 
yields, as discussed above. An old er grid of yields is also plot- 
ted, which consisted of the WW95 SNII yields up to 25 Mq and 
the SNIb/c yields associated with the MM97. non-rotating WR 
models. For these SNIb/c yields, the yield versus CO core mass 
relation was derived from the pure He star models o f lWLW95L 
Regardi ng the -''Al SN yields of WH07, they are very similar to 
those of lLC06l up to 40 (variations remain within ± 50%) ex- 
cept for the stars around 14-15 M©, for which ,LC06. fou nd a dra- 
matic increase of th e ^^Al S N yield not seen i n IWH07l . Beyond 
40 M0, the yields of lWH07l are below those of lLC06l bv a factor 
of up to 10 for the I2OM0 model. This is very like ly an effect 
of the stronger mass loss experienced by the lWHOTl higher mass 
models. 

Again, n o detailed comparison w as performed but we 
noted that the lWellstein & Langed (119991) mass lossprescription 
(scaled by 1/3 to account for clumping) used by |WH07 gives 
mass loss rates in the late WR stages (t ypically WNE and WCO 
stars) that are a factor of 2-3 above the lNugis & LamersI (l2000h 
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Fig. 4. <'"Fe supernova yields fro m [LCOel and IWHOTl Also plot- 
ted are th e yield s for IMMP05I extrapolated through CO core 
mass from lLCQg 



prescription used by lLC06l This is obyious from the final presu- 
pernova masses of the 40-120 models: while those obtained 
by LC06 are in the 12-2OM0 range, those obtained by WH07 
are between 6-IOM0. It is interesting to note that LC06 tested 
the impact of increased mass loss during t he WN E and WCO 
phases by using the prescription of Lanaej ( 19891). which giye s 
rates a factor of 2-3 above the rates of Nugis & Lamers ("2000). 
In this case, their 40-1 20 Mq models ended up with final masses 
and ^^Al SN yields that are very similai" to those of IWHOTi A 
stronger mass loss actually reduces the size of the He and CO 
cor es and l eads to lower ^*A1 SN yields. The difference between 
the IWHOTi and LC06 ^^Al SN yields is all the more important 
because in the WH07 yields there is a contribution from the v- 
process tha t is abse nt in the|LC06 yields (see l4.1.3l l. The alterna- 
tive grid oflWHOTi thus illustrates how important mass loss can 
be for nucleosynthesis. 

In Fig.H the ''"Fe yields of LC06 and WH07 are plotted with 
the SN yields extrapolated through the CO core mass for the 
ImM POS rotating WR models. As for ^""Al SN yields, the eff'ect 
of rotation on the size of the CO core, hence on the extrapolated 
yields, is manifest. The WH07 yields are apparently quite dif- 
ferent from those of|LC06j. Between 10-20 they are up to a 
factor of 20 above tho se of LC06, then between 20-50 they 
are still above those of iLCOd bv more moderate factors of up to 
5, and beyond 50 M© they drop below LC06 by quite large fac- 
tors. For the higher mass models, the lower yields of WH07 very 
likely come from the stronger mass loss of WR stars, as a lready 
discussed above. Here again, the same effect was noted bvlLC06l 
when they used the increased mass loss rates of lLangeJ (Il989l) . 
For the lower masses, however, the differences are not so easy 
to explain. The lower initial metallicity should have resulted, all 
other parameters being unchanged, in lower yields since ^°Fe is a 
pure secondary product. The major reaction rates for the produc- 
tion and destruction of ^"Fe were found to be identical in WHOj 
and ILC0 6, so the difference may come from another aspect of 
stellar physics such as convection or mixing. 



3.3. The stellar content of the Cygnus complex 

The last ingredient of our simulation is the definition of the 
stellar population to be modelled. The peculiarities of the 



Cygnus region have fostered a great number of studies that 
together provide us with a pretty good knowledge of the area. In 
particular, recent works have lead to a substantial revision of the 
stellar content of the Cyg OB2 cluster and thus allow for more 
accurate predictions. 

From the stellar census of iKnodlsedei 



iKnod lseder et alJ (l2002l) . and 'Le Duigou & Knodlseder' (2002 
we selected the most relevant OB associations and open clusters 
of the Cygnus region. We restricted ourselves to a longitude 
range of 70-90° and among the stellar groups falling in that 
interval, we excluded the ones with ages above 20Myr (since 
their decay fluxes are negligible after that time, see Fig. |6]l. 
The selected objects are listed in Table [T] t ogether with their 
main characteristics. IKnodlseder et alj (l2002h indicate that most 
open clusters are potentially physically related to some of the 
OB associations. In total, our selection amounts to a total of ~ 
170 O stars, most of which are in the Cyg OB2 cluster. This 
ensemble is what we have called so far the Cygnus complex. 
ISchneider et al.l (l2006l) demonstrate that Cyg OB 1 , 2, and 9 very 
likely originated from the same giant molecular cloud, which 
gives a physical reality to this grouping. The connection of Cyg 
OB3, 7, and 8 to this complex is, however, still uncertain. 

Whether all these clusters share a common origin does not 
matter at all for our purpose. Our aim here is simply to gather the 
nearby clustered stellar population that is very likely responsible 
for most of the stron g 1809keV emis s ion from Cygnus. In s ome 
previous studies (fdeT Rio et al. I ll996t iPliischke et al.ll200l . the 
total 1809keV flux from the Cygnus region was compared to 
theoretical predictions based on the clustered stellar content of 
the Cygnus region supplemented by a few isolated WR stars 
and SNRs. In our study, we focus both observationally and 
theoretically on the clustered stellar content alone, because the 
clustered population is likely to be bette r invent oried than the 
isolated or non-clustered population. In iPaper iL we used the 
INTEGRAL/SPI data to disentangle the 1809keV flux of the 
Cygnus complex from the more spatially extended emission 
likely due to unclustered objects like non-member WR and 
OB stars (runaway stars or stars from previous star formation 
episodes, see Comeron & Pasquali (2007) and IComeron et al.l 
(1200 8)) or isolated SNRs. In the present paper, we compare 
this flux to population synthesis predictions based only on the 
Cygnus complex population. 

3.4. The code 

Our population sy nthesis code is very similar to many other such 
tools (Uke Leither eretaT]ll992t ICervino et ai]l20od IVoss et al.1 
2009, for example). For the case we are interested in, the param- 
eters of a simulation are the characteristics of the stellar popula- 
tion to be modelled: slope and bounds of the IMF (initial mass 
function, assumed here to be a homogeneous power-law func- 
tion), stellar content (given as the number of stars observed in a 
certain mass interval), and mean distance to the cluster (to com- 
pute the decay fluxes). We assume that the star formation occurs 
as an instantaneous burst, which means that all stars form at the 
same time (discussed later). 

Once the population to be synthesised has been defined, we 
create a cluster realisation by randomly sampling the IMF un- 
til the observational criterion (the number of observed stars) is 
met. Then, all stars are evolved in time steps typically of 10"* yrs. 
For each run, we compute the time-dependent release of ^^Al 
and "'Fe, the decay fluxes at 1809 and 1 173/1332keV and a few 
other outputs such as the mechanical luminosity. To evaluate the 
effects of finite sampling, a given stellar population is modelled 
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Table 1. Characteristics of the Cygnus OB associations and stel- 
lar clusters collectively refeiTed to as the Cygnus complex (from 
[Knodlseder et al. 2002). 



Name 


Observed population 


Distance (pc) 


Age (Myr) 


Cyg OBI 


23 6 [15,40 Mq] 


1905 


4 


Cyg OB2 


120 e [20,120 Mq] 


1584 


2.5 


Cyg OB 3 


14 6 [25,60 Mo] 


2187 


3.5 


Cyg OB7 


10 € [7,25 Mq] 


832 


3.5 


Cyg OB 8 


6 s [20,40 Mq] 


2399 


7.5 


Cyg OB9 


8 6 [20,40 Mo] 


1259 


3.5 


Ber86 


11 6 [7,25 Mo] 


1660 


4 


Ber87 


24 6 [7,25 Mo] 


1905 


4.5 


NGC 6871 


13 6 [7,25 Mo] 


2399 


5.5 


NGC 6913 


13 6 [7,40 Mo] 


1820 


3.5 


NGC 6910 


7 6 [7,25 Mo] 


1820 


4.5 


NGC 6883 


2 6 [12,15 Mo] 


1820 


15 


IC 4996 


6 6 [6,25 Mq] 


1660 


5.5 



several times (typically 100-200) until the mean stellar mass and 
number of stars obtained over all trials are within 1 % of the the- 
oretical values for the specified IMF. For each quantity, this al- 
lows computing the mean value and the statistical variance at 
each time step in the life of the cluster. 

Although a synthetic star can take any initial mass between 
11 and I2OM0 (the bounds of our grid of yields), we have the 
nucleosynthesis yields for only a finite number of initial stel- 
lar masses, t ypically 1 1, 12, 13, 14, 15, 16, 17, 20, 40, 60, 80, 
120 Mo (the IMMP05I grids of the higher mass models differ 
for the various metallicities). We therefore need an interpola- 
tion scheme to compute the yields of the synthetic stars from the 
ones at our disposal. 

The characteristics of a synthetic star with a random initial 
mass are determined from the closest lower and higher mass 
models (for instance, the characteristics of a 53.6 Mq star are 
determined from the 40 and 60 M© models). A first cubic-spline 
interpolation is performed to determine the durations of the main 
evolutionary sequences of the star: central H-burning, central 
He-burning and final stages (when available). Then, for each of 
these three sequences, the ^''Al and ''"Fe quantities ejected over a 
given time fraction of the sequence are obtained by linear inter- 
polation of the quantities ejected by the closest lower and higher 
mass models over the same time fraction. When the synthetic 
star explodes, its ^^Al and ^°Fe yields are determined by a linear 
interpolation based on the CO core mass and released in the cur- 
rent time step. 

The work of lLC06l clearly shows that the yields are not 
smooth functions of the initial stellar mas^. For instance, the 
lack of a C convective shell in their 14 and 15 M© models trans- 
lates into a dramatic increase in the Al production compared to 

On the contrary, the durations of the main burning stages are a quite 
smooth function of the initial mass, close to a power-law, and that is why 
in our work they are determined by a cubic-spline interpolation instead 
of a simple linear interpolation. 



the 13 or 16 M© models. As a consequence, it should be remem- 
bered that the finite grid of yields and the associated interpola- 
tion scheme are potential sources of errors. 



4. Predictions versus observations 

Based on the code and inputs presented in Sect. [3] we have sim- 
ulated the history of the ejection of ^''Al and ^°Fe in the Cygnus 
complex and computed the corresponding decay fluxes. We first 
focused on the purely photometric aspect, trying to reproduce 
the observed fluxes by simply adding the contributions from all 
the stellar groups listed in Table [T] In the second part, we also 
attempted to reproduce the spatial extent of this emission. 

4. 1 . The decay emission fluxes 

All the stellar groups listed in Table [T] have been modelled in- 
dependently, assuming a coeval star formation for each of them. 
Then, the lightcurves at 1809 and 1173/1332keV were added, 
with the appropriate time shifts to take the differences in age 
of the different stellar populations into account. The resulting 
lightcurves for the whole Cygnus complex are given with t=0 
corresponding to the present day and thus allowed exploring the 
past and future evolution of the decay fluxes. 

As a first step, we focused on the results obtained for an ini- 
tial solar metallicity, firstly because most of the past population 
synthesis studies were based on nucleosynthesis yields from so- 
lar metallicity stellar models so a comparison could be made, but 
also because the Cygnus complex lies on the solar circle and is 
thus expected to have a nearly solar metallicity. 

4.1 .1 . The case for solar metallicity Z=0.02 

The results obtained for an initial solar metallicity are pre- 
sented in Fig.|5](top), together with the observational constraints 
we derived from the INTEGRAL/SPI observations. The pre- 
dicted decay fluxes at 1809 and 1173/1332keV are there to- 
gether with their typical variances due to the finite IMF sam- 
pling. Fro m the open clusters and OB associations ages deter- 
mined by iKnodlseder et al.l (l2002h . our current position along 
the time axis is indicated by the zero point. Yet, these clus- 
ter ages were determined by fitting isochrones to Hertzprung- 
Russel diagrams, and the isochrones used by these authors were 
built from the non-rotating stellar models available at that time. 
Rotation is now known to increase the lifetime of the models by 
15-2 5%, mainly through augmenting the main-sequence dura- 
tion (Mevnet & Maeder 2003). Therefore, the cluster ages have 
probably been underestimated by the same amount, and our cur- 
rent position on the time axis of Fig.|5]is likely shifted to the 
right of the zero-point, as indicated by the grey shaded area. 

Inside this time range, there is clear agreement between the 
predicted 1809keV flux from ^^Al and the observed value, while 
the theoretical level of the ^"Fe decay emission is consistent 
with our upper limit. Our population synthesis indicates that we 
are currently in a period of steep increase of the 1809keV flux 
(over Myr timescales) that follows from a substantial release of 
^^Al through the stellar winds of the most massive stars, sup- 
plemented by the first supernova explosions that occurred in the 
oldest stellar groups, typically those above 3.5 Myr (the short- 
est lifetime in our grid of models, corresponding to an initial 
120 Mq star). At about 1 Myr on the time axis of Fig. |5](which 
may be our current position), most of the wind ^^Al has been 
injected in the ISM; the supernovae then take over and now re- 
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lease both ^^Al and *°Fe. From there on, the 1809keV flux is 
expected to rise again, less rapidly though, and to reach within 
~ 1 Myr a maximum mean value of ~ 4.2 x lO^^phcm^^s"', 
which remains in our measured flux interval. Regarding the 
1173/1332keV emission, the most massive stars of the oldest 
stellar groups exploded and consequently started injecting *"Fe 
in the ISM. The number of events up to now, however, is too low 
to have built up a detectable amount of *"Fe. Moreover, our pre- 
diction shows that even the maximum 1 173/1332keV flux, to be 
reached in ~ 2-3 Myr from now, remains below our present up- 
per limit. 

Our simulation predicts that a certain number of SNe, typ- 
ically 10-20 over the last Myr, exploded in the Cygnus com- 
plex, which contradicts the apparent absence of radio SNRs and 
the detection of only a single pulsar in its various OB associ- 
ations and open clusters. The absence of SNRs could be ex- 
plained by the fact that SNRs expanding inside the hot tenuous 
interiors of superbubbl es exhibit v ery weak characteristic radio 
and optical signatures (IChulll99"7h . With their high stellar con- 
tent, the OB associations of the Cygnus complex are expected 
to blow large superbubbles inside of which the SNR are likely 
to be invisible. Moreover, the typical radio lifetime of a SNR 
is short compar ed to the lifetime of ^^Al (10** versus 10''yrs, 
iFrail et al.|[T994b and all SNRs may well be completely diluted 
in the ISM by now. In that case, we would expect to see some 
manifestations of their compact remnants. Up to now, only a sin- 
gle pulsar has been proposed associat ed with the Cygnu s com- 
plex, more precisely with Cyg OB2 dCamilo et al.l l2009i) . de- 
spite the dedicate d deep surveys that have been undertaken at 
radio frequencies (iJanssen et alj|2009l) . The difficulty of finding 
pulsars in the Cygnus complex could come from an excessive 
scintillation or dispersion of the radio signals in the strongly 
ionised environment of such a concentration of massive stars. 
Or it may be that the narrow cones of their radio emission do 
not intersect our line of sight. In the latter case, these pulsars, 
if they exist, may become "visible" to the Fermi/LAT instru- 
ment through their gamma-ray emission, which is unabsorbed 
and less-focused compared to the radio emission. The discov- 
ery by Fermi/LAT of the gamma-ray pulsar PSR J2032-I-4127 
that is li kely located in Cyg OB2 is promising in this respect 
dAbdo et al. 2009) . 

The above result constitutes a step forward compared to the 
previous studies of the nucleosynthesis activity of the Cygnus re- 
gion. Indeed, the latest evaluations of the situation conclude that, 
from the established stellar content (including the revised 100- 
120 O stars of Cyg OB 2) and the stell ar yields of iMevnet et alJ 
(Il997h and IWooslev & Weaved (Il995h . the observed 1809keV 
flux from Cygnus is a factor of 2-3 above the predictions (see 
for instance Knodlseder et al. 2002). Pluschke et al. (2000) point 
out that the possible enhanced yields of massive close binary 
systems could alleviate the discrepancy. While it is clear that 
such systems may occur in the Cygnus complex and could in- 
deed have a strong impact on the nucleosynthesis yields, our 
theoretical knowledge of this scenario is still too uncertain to 
allow us to take it quantitatively into consideration. Following 
the significant upward revision of the Cyg OB2 stellar con- 
tent ( Knodlseder 2000). lPluschke et all (120021) suggested that the 
stellar content of the other OB associations of the Cygnus re- 
gion may be underestimated due to a strong interstellar extinc- 
tion; with correction factors based on the CO column densities 
in the direction of Cygnus, the authors show that the resulting 
increase in the stellar content could solve the problem, and yet, 
this revision of the stellar population ha s not been observation- 
ally substantiated since then. Moreover. ISchneider et al.l (l2006l) 



show that some fraction of the CO mass observed toward Cygnus 
is actually located behind the OB associations, which means that 
the applied correction was incorrect in some cases. 

Our work indicates that the past dilemma resulted from an 
overestimate of the 1809keV flux of the Cygnus complex, com- 
bined with an und erestimate of the nucleosynthesis yields. First 
of all, in Paper I, a careful analysis of the INTEGRAL/SPI 
1809keV observations showed that the decay flux attributable 
to the Cygnus complex is only ~ 65% of the total flux coming 
from the Cygnus region. This result was obtained by separating 
the emission due to the Cygnus complex from the foreground 
and background mean Galactic contribution. Then, stellar rota- 
tion, as discussed in Sect. 13.11 turned out to be a key factor in 
nucleosynthesis since it enhances the production of Al but also 
affects the rate of its release in the ISM. In addition, the compu- 
tation in a coherent way of the SNIb/c contribution (by LCO^ 
showed that the SN yields of the most massive stars cannot be 
neglected or simply extrapolated from pure He stars calculations. 

4.1 .2. The case for metallicity Z=0.01 

Most of the past studies were based on the nucleosynthesis yields 
of solar metallicity stellar models, so we presented the above re- 
sult for comparison; yet, there are indications that the stellar pop- 
ulation of the Cygnus region may be subsolar with a metallicity 
of about Z=0.01 (Daflon et al,. 20Qi ). Although we do not have a 
fully homogeneous grid of stellar models of various metallicities 
covering hydrostatic and explosive phases at our disposal, we ex- 
ploited as much as possible the works of LC06 and MMP05 to 
estimate the ^^Al and ''"Fe yields of a non-solar population of 
massive stars. We recall here that our grid of yields was built 
from purely hydrostatic stellar tracks of various initial metallic- 
ities combined with full calculations (hydrostatic and explosive) 
only at solar metallicity, the connection being made (when nec- 
essary) through the CO core mass (see l3.2l ). 

The decay lightcurves for Z=0.01 (obtained by linear inter- 
polation between Z=0.008 and Z=0.02) are presented in Fig. 
II] (bottom). Compared to the Z=0.02 case, the peak 1809 keV 
flux is somewhat reduced and occurs at a later time, while the 
1173/1332keV flux is increased by almost a factor of 3. Both 
differences come from the effect of lower metallicity on stel- 
lar evolution. Indeed, a lower initial metallicity implies reduced 
mass loss for stars initially more massive than 20 Mq, since the 
winds of the OB and WR phases are mostly driven by radia- 
tive pressure on metals through their numerous absorption lines. 
Reduced mass loss leads to larger CO core masses, hence higher 
SN yields, as we assumed that the latter mostly depends on the 
CO core mass. This accounts for the considerable increase in the 
''"Fe production in the Z=0.01 case. For ^^Al, however, the in- 
crease in the SN yields is compensated by a decrease in the wind 
yields; indeed, a lower initial metallicity means a reduced mass 
loss but also a reduced initial ^^Mg reservoir for proton cap- 
ture (hence synthesis of ^^Al) during central H-burning. On the 
whole, the 1809 keV lightcurve is not dramatically altered by the 
reduced metallicity, but the stronger SN contribution explains 
the delay in the 1809keV peak flux compared to the Z=0.02 
case. 

Within the likely present-day time window, the mean pre- 
dicted 1809keV flux is at best a factor of two below the ob- 
served value. Consistency is found later on, at ~ 2-3 Myr from 
now, but then the 1 173/1332keV flux has risen close to or even 
above our upper limit. It therefore seems that the overall agree- 
ment obtained at Z=0.02 does no longer holds if the subsolar 
metallicity of Cygnus is accounted for We may question this 
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subsolar metallicity. It is indeed somewhat surprising that stel- 
lar associations lying at about the same galactocentric distance 
as the Sun and be ing much younger could be subsolar, although 
iDaflon et alJ (12001 ) mention that their abundances are consis- 
tent with the expected values, as predicted by various Galactic 
abundance gradients. While the C, N O, Si, Al, and Fe were 
found to be underabundant relative to solar values, Mg and S 
were found to be solar and Mg has straightforward relevance for 
^^Al nucleosynthesis. The authors state that B stars in the solar 
neighbourhood generally have subsolar abundances, so one may 
worry about a general systematic bias in determining chemical 
abundance s in B stars. We also note that a solar abundance was 
derived by iNajarrol (1200 ih for the P-Cygni star, which lies at ~ 
2 kpc, close to the Cygnus complex. 

From these arguments, we therefore consider that the metal- 
licity of the Cygnus complex is s till uncertain , especially since 
the associations investigated by iDaflon et alj (1200 1) are Cyg 
OB3 and 7, which are not those contributing the most to the 
Cygnus complex. In addition to this, we discuss below several 
other potential sources of error or uncertainty that can mod- 
ify and/or alleviate the current situation for the gamma-ray line 
emission of the Cygnus complex. 

4.1.3. Sources of error/uncertainty 

Up to now, we have only considered the uncertainty on the pre- 
dicted fluxes due to finite IMF sampling, but the uncertainty 
on the exact distances and ages of the stellar associations of 
the Cygnu s complex should also be taken into account. From 
iKnodlseder et al.l (120021) . we estimated that around the present 
date, the uncertainties on the ages and distances of the stellar 
clusters translate into typical uncertainties of ~ 1.0 and 0.25 x 
10"^ phcm"^ s"' on the decay line fluxes of ^^Al and ^'^Fe, re- 
spectively. In connection with the age of the stellar clusters, it 
should be noted that our assumption of a coeval star formation 
favours maximum peak fluxes. The most massive stars (between 
30 and I2OM0) have very similar lifetimes, so in a coeval popu- 
lation the nucleosynthesis contributions of their last stages (WR 
and SNe) are released at almost the same time and therefore lead 
to a steep rise in the decay fluxes. In contrast, if star formation 
occurs over a few Myrs, the nucleosynthesis contributions are 
spread over a longer time interval and the resulting peak fluxes 
are lower. 

We emphasise once again that our grid of nucleosynthesis 
yields is not homogeneous and self-consistent and that the SN 
yields of the higher-mass stars at all metallicities were extrap- 
olated from non-rotating solar models (whereas the hydrostatic 
phases of these higher mass stars were computed by including 
stellar rotation). As a first potential issue, the direct effect of 
metallicity on the SN yields of non-solar models is not taken into 
account. Indeed, our grid of yields takes into account the effect of 
metallicity on stellar evolution processes like the mass loss and 
the late stage core sizes. But the direct influence of metallicity, 
that is the differing chemical composition of the stellar matter 
from which the nucleosynthesis proceeds, is not considered. For 
instance, we have shown before that an initial subsolar metal- 
licity leads to bigger CO core, hence presumably to higher ''"Fe 
SN yields; however, this should be mitigated by the fact that a 
lower metallicity also means reduced initial quantities of Fe and 
CNO for the synthesis of ^"Fe. Since ''"Fe is a pure secondary el- 
ement, it is likely that our grid overestimates the ''"Fe SN yields 
of the subsolar metallicity models. Regarding ^''Al, the direct 
impact of the metallicity is less straightforward. Indeed, most of 
the ^''Al released by supernovae is of explosive origin and the 
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Fig. 5. Predicted versus observed decay fluxes, for Z=0.02 (so- 
lar) initial metallicity (top) and Z=0.01 initial metallicity (bot- 
tom). The red curve corresponds to the 1809keV flux and the 
blue curve to the 1173/1332keV flux. The dotted lines give the 
typical variance due to the finite IMF sampling. The orange band 
indicates the observed 1 809 ke V flux and the green line the upper 
limit on the 1 173/1332keV emission from the Cygnus complex. 
The vertical shaded area indicates our probable current position 
along the time axis (see text). 



corresponding nuclear processes seem to depend only mildly on 
the initial metallicity (although the many involved species do not 
allow drawing definitive conclusions). In this case, our original 
assumption that the SN yields predominantly depend on the CO 
core mass appears to be more reliable for ^^Al than for ''"Fe. 

A second potential issue stems from SN yields coming from 
non-rotating models. Rotation has proven to be efficient at en- 
hancing the ^^Al production during central H-burning, mostly 
thanks to an increased convection and mixing. If rotation were to 
modify the yields of the C/Ne shell-burning episodes in a similar 
manner, the corresponding ^^ Al and ^"Fe production would be 
higher than those computed bv lLC06l In addition, increasing the 
shell convection efficiency would also alter the final mass-radius 
relation on which our estimate of the explosive nucleosynthe- 
sis depends. We should also mention that, in the grid of yields 
we are using, the v-proces^ has not been taken into account. 



' In the v-process, protons liberated by the v-spallation of e.g. -"Ne, 
"'O, ^""Na, and -''Mg capture on -'Mg to form ^*'A1. 
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This production channe l was estimated to enhance the explosive 
^'^Al y ields by 30-50% (IWooslev & Weav eijll995l;lTimmes et all 
These values may have to be revised with more sophisti- 
cated and realistic neutrino spectra and transport. 

More generally, the nucleosynthesis yields should be re- 
garded as sensitive to most ingredients of stellar models, such 
as convection and mixing, opacities or nuclear reaction rates. In 
13.21 we compared the LC06 and WH07 sets of yields and showed 
that the choice of mass loss prescription or the implementation of 
convective physics does significantly affect the nucleosynthesis 
of Al and *"Fe. We did not observe drastic variations by many 
orders of magnitude in the outputs of the two grids of stellar 
models, but the typical differences are high enough to possibly 
modify our conclusions. Therefore, although they were not in- 
cluded in our population synthesis, we qualitatively assessed the 
impact that the WH07 yields would have on the predicted decay 
fluxes from the Cygnus complex. 

For stars in the 20-1 20 Mq mass range (which dominate the 
current nucleosynthesis activit y of th e Cygnus complex), the 
^^Al wind and SN yields of the'WHOT" models are clearly lower 
than those of MMP05 and LC06 (see Figs. |2] and [3]). We esti- 
mated that the use of these yields as such in our population syn- 
thesis would lower the 1 809 ke V flux by at least a factor of 2 and 
thus lead to a mismatch between predictions and observations if 
the Cygnus complex has nearly solar metallicity. Nevertheless, 
the WH07 models do not include stellar rotation, which turned 
out to strongly enhance the ^^ Al yields . Then, keeping the hope- 
fully more realistic yields of |MMP()5I fo r the ^^ Al wind contri- 
bution, the lower ^*A1 SN yields of the |WH07| models are not 
expected to seriously alter the rising part of the early 1809keV 
lightcurve that we observe today (see Fig.|5]l, because the latter 
is dominated by the wind contribution. This is clearly apparent 
in Fig.|6]for the simulation of Cyg OB2 only: when the first stars 
of the coeval cluster explode, at ~3.5Myr, the flux has already 
risen to 70-80% of its maximum value and is already consistent 
with the observed value. Consequently, even if the WH07 mod- 
els may lead to less ^*A1 being ejected by SNe, the wind yields 
of the presently available rotating models alone can account for 
most of the observed ^*A1 decay flux. Regarding ^''Fe, the WH03 
yields are not very diff'erent from the yields we have been us- 
ing (MMP05 models extrapolated through CO core mass from 
iLCOfifi " over the 20-120 Mq mass range (see Fig.|4ll. We thus ex- 
pect the 1 173/1332keV lightcurve to still be consistent with our 
upper limit event if the WH07 yields were used. 

We therefore believe that the WH07 set of yields would not 
modify our results for the solar metallicity case if we still use the 
enhanced ^^Al wind yields of the rotating m odels. A definitive 
conclusion would require fully implementing'WH07 in our pop- 
ulation synthesis code (in order, for instance, to take the specific 
SN yields to CO core mass relation of this grid properly into 
account). More generally, as alternative sets of yields become 
available, especially from rotating models, it will be interesting 
to see if our conclusions still hold. Even in the absence of ma- 
jor improvements of the stellar models, the evolution of some of 
their inputs alone, such as nuclear data, can alter the nucleosyn- 
thesis results. Some critical reaction rates, such as a(2a,jy^C 
or '^C(Q',y)"'0, are k nown to have strong effects on the final 
yields (see for instance FTur et al.P 2009). As an illustration of the 
uncertainties on nuclear data, the half-life of *"Fe was recently 
revised to 2.6 Myr, instead of the 1 .5 Myr val ue we have been 
using throughout this work (iRugel et alJl2009h . In terms of pre- 
dicted 1 173/1332keV lightcurves for the Cygnus complex, how- 
ever, this longer lifetime simply implies an even lower emission 
at early times, still consistent with our upper limit (for a more 



quantitative ana lysis of the impact of the new ^°Fe lifetime, see 
IVoss et al.ll2009 '). 

We also emphasise that all the result s presented so far are 
based on a "classical" IMF slope of 1.35 (ISalpeteilll955h . Since 
the universality of this parameter is still a matter of debate, we 
show in Fig.|6]how the decay lightcurves for Cyg OB2 vary with 
the slope of the IMF. The impact is fairly straightforward: a flat- 
ter IMF means more very massive stars (for a given total number 
of massive stars) and the 1809keV peak at 3-4 Myr attributable 
to the winds and SNe of WR stars is consequently stronger, while 
a steeper IMF means fewer very massive stars, which reduces the 
1809 peak flux at 3-4 Myr and enhances the flux at 15-16Myr, as 
a result of more SN explosions of 13-15 M0 stars. 

In connection with the potential effect of the IMF slope, we 
want to draw attention to the possible impact of binarity. It is 
commonly accepted that for stellar clusters more distant than a 
few kpc, most of the observed point-sources assumed to be mem- 
bers of the clusters may actually be unresolved binaries or mul- 
tiple systems. The census of their stellar content is consequently 
biased, particularly for purely photometric approaches, but even 
when relying on spectroscopy. Kobulnickv & Frver (2007) find 
that the binary fraction of Cyg OB2 is very high, and may well 
be close to 1 ; moreover, it seems that massive stars preferentially 
have massive companions. We therefore wonder about the pos- 
sibility that a substantial fraction of the 100-120 O stars of Cyg 
OB2, many of which have been identified photometrically from 
the 2MASS data, may in fact be unresolved massive binaries or 
multiple systems. If this were true, the number of massive stars 
in our population synthesis would have to be increased, which 
could lead to higher predicted fluxes. To properly account for 
this binarity effect, however, the various observations devoted to 
Cyg OB2 (and the other OB associations) have to be thoroughly 
compared to clearly identify which fraction of the objects as- 
sumed today as single massive stars may actually hide multiples, 
but this is beyond the scope of the present paper. 



4.2. The spatial diffusion of^^AI 

In the previous sections, we made use of the photometric aspects 
of the INTEGRAL/SPI observations of Cygnus alone. Here we 
want to extend our interpretation of the SPI results to the spatial 
and spectrometric information provided by this instrument. 

The numerous massive stars of the Cygnus complex are 
expected to have blown one or several superbubbles (SBs) 
through the combin e d action of their win d s and SN explosions 
(Castor et al.' l9T?; 



19771: 'Tomisaka & Ikeuch 



Weaver et al 
1986; Mac Low & McCray 1988[ 
Soft X-ray observations by HEAO-3 and ROSAT have revealed 
a large structure in the Cygnus direction that may well be the 
outermost part of a gigantic SB, the centre of \yhich is concealed 
by intervening interstellar material dCash et al.ll9 80'). Fresh ^^Al 
is expected to be injected inside the hot and tenuous interior of 
this SB, and its diffusion on large scales will then be driven by 
the growth of that structure. 

Under certain simplifications, a self-si milar analytical sol u- 
tion can be found for the growth of the SB dCastor et al.lll975h : 



l/5„-l/5,3/5 



'38 



pc 



■ 40xLli'n;'''tf'' kms-' 



(1) 
(2) 



where Rs and Vs are the outer radius and velocity of the SB, 
respectively, given as functions of the mechanical luminosity of 
the central cluster L38 (in 10^^ erg s" ' ), the uniform density of the 
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Fig. 6. Decay fluxes of Al and Fe from Cyg OB2 as a func- 
tion of the IMF slope for an initial solar metallicity. In contrast 
to Fig.|5] the t=0 point does not mark the present day but instead 
the date of the coeval star formation. 



surrounding medium «i (in cm"^) and the age of the central clus- 
ter t(, (in Myr). These equations will allow first estimates of the 
expected size of the 1809keV diffuse emission from Cygnus, but 
also show the relatively weak dependence of the bubble's size on 
some uncertain parameters, such as the density of the surround- 
ing ISM medium. 

In the following, we focus on the Cyg OB2 cluster only 
because we expect most of the mechanical power and nucle- 
osynthesis activity of the Cygnus complex to come from it. 
From our population synthesis tool, we have found that the 
mechanical luminosity of Cyg OB2 is about 4 x 10^** ergs"' 
over the first 3 Myr of the cluster's lif^ Using this value, to- 
gether with a typical ISM density of 1 cm we find that the 
radius of the SB blown by Cyg OB2 is about 170pc after 3 Myr 
Assuming that the ^^Al homogeneously fills the interior of the 
SB, the corresponding 1809keV decay signal would appear to 
us as a diff'use emission with a typical size of 12° (given the 
~1.6kpc distance of Cyg OB2). This estimate tallies with our 



* Mor e details ab out the calculation of the mechanical power can be 
found in lVoss et al.l (2009), who also developed a population synthesis 
code dedicated to gamma-ray line flux predictions. Their code is based 
on similar assumptions and the same set of stellar data, and the results 
of both tools have been checked to be consistent. 



INTEGRAL/SPI measurement of a 9-10° angular size (espe- 
cially since the maximum size of the emission is weakly con- 
strained by SPI, see l2.1l i. It is also consistent with previ ous mea- 
surements by CGRO/COMPTEL (jPliischke et allfeoOlh . 

We advanced the accuracy of our estimate of the expected 
extent of the 1809keV diffuse emission from Cygnus by using a 
numerical simulation. Our observational constraints are actually 
not tight enough to justify such an effort, but the present study 
was an opportunity to address the question of the diffusion of 
^^Al around a stellar cluster in some details. We used the VH- 
I hydrodynamical cod^ to compute the evolution of a simple 
spherical bubble blown by the injection of kinetic energy at the 
centre. The code was modified to include an equation for the ad- 
vection and decay of Al in addition to the three classical Euler 
equations. A variable input was also implemented to allow injec- 
tion of -^Al according to the time profile given by our popula- 
tion synthesis code. Last, radiative losses, which rapidly become 
significant in the swept-up ISM shell, were implemented in the 
form of a cooling function from I Sutherland & Dopital (Il993h . 

By again using a constant mechanical luminosity of 4 x 
10^^ ergs"' with a typical wind velocity of 2000 km s"', and an 
ISM density of 1 cm"^ for a medium with 90% H and 10% He (in 
number), the resulting ^^Al mass distribution after about 3 Myr 
is plotted in Fig. Q The inner flat distribution corresponds to 
the rapid advection with the freely-expanding winds, followed 
by a sharp discontinuity that marks the inner boundary of the 
shocked-wind region that actually constitutes most of the inte- 
rior of the SB. A sudden drop in the ^^Al content then occurs at 
the contact discontinuity that separates the shocked stellar mate- 
rial from the shocked swept-up ISM. From Fig.|7] it appears that 
the ^''Al fills the entire SB but somehow accumulates at roughly 
half its radius. The ^^Al mass distribution is relatively homoge- 
neous inside the hot shocked interior of the SB, with variations 
along the radius that remain within a factor of 10. 

According to our numerical simulation, the angular size of 
the 1809keV emission from the -^Al released by the Cyg OB2 
cluster is ~ 10°, which is again consistent with the measure- 
ments by INTEGRAL/SPI and CGRO/COMPTEL. Our simula- 
tion also indicates that most of the -*A1 is aff'ected by expansion 
velocities in the 50-300 km s"' range, which agrees with the ve- 
locities inferred from the 1809keV line broadening. We want to 
caution, however, that several physical processes were neglected 
in our numerical model. First of all, thermal conduction between 
the hot tenuous interior of the SB and the cold dense swept- 
up ISM shell was not taken into account (ICastor et al.l 119751: 
I Weaver et al.lll97'7b . This energy transfer causes material from 
the shell to evaporate and enrich the mass content of the SB (the 
resulting mass flow actually dominates the mass carried by the 
stellar winds). Nevertheless, Tomisaka & Ikeuchi ( 1986), show 
that this process does not significantly alter the overall size of 
the structure, but rather acts as an internal mass-energy redistri- 
bution. Second, our one-dimensional simulation cannot account 
for the turbulence that very likely occurs in the SB as a result 
of the multiple wind flows inside the stellar cluster, inhomo- 
geneities in the surrounding medium or hydrodynamical insta- 
bilities. Such turbulence will mix the stellar material inside the 
bubble and probably ablate mass from the swept-up ISM shell. 
These two effects will obviously blur the contact discontinuity 
and alter the mass distribution compared to what we presented 
in Fig.|7] Turbulent motions will also be an additional source for 
the broadening of the decay line. 



' Made available to the community by John Blondin and collabora- 
tors, see http://wonka.physics.ncsu.edu/pub/VH-l I 
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Fig. 7. Distribution of Al inside the superbubble blown by Cyg 
OB2 after 3 Myr, as predicted from a simple numerical simula- 
tion. 



5. Conclusion 

We have compared the INTEGRAL/SPI observations of the 
gamma-ray line emission of the young OB associations and open 
clusters of the Cygnus region with theoretical expectations based 
on the latest models of stellar nucleosynthesis. This nearby com- 
plex of massive stars exhibits an unambiguous diffuse emission 
at 1809keV caused by the decay of ^^Al but no detectable sig- 
nal of the decay of ^"Fe at 1 173/1332keV. The latest studies of 
the nucleosynthesis activity of the Cygnus region show that the 
1809keV emission deduced from CGRO/COMPTEL observa- 
tions was a factor of 2-3 above the predictions from the theoret- 
ical yields available then and the observed stellar content of the 
Cygnus region. 

Compared to previous works, our grid of yields includes 
some of the effects of stellar rotation for the higher-mass stars 
and a coherent estimate of the contribution from SNlb/c. From 
this, we show that the observed decay fluxes from the Cygnus 
complex are consistent with the values predicted by popula- 
tion synthesis at solar metallicity. Our work indicates that the 
past dilemma resulted from overestimating the 1809keV flux of 
the Cygnus complex, combined with underestimating the nucle- 
osynthesis yields. We showed in a previous work, through a care- 
ful analysis of the INTEGRAL/SPl 1809keV observations, that 
the decay flux attributable to the Cygnus complex is only ~ 65% 
of the total flux coming from the Cygnus region. This result was 
obtained by separating the emission due to the Cygnus complex 
from the foreground and background mean Galactic contribu- 
tion. Then, stellar rotation turned out to be a key factor in nu- 
cleosynthesis since it enhances the production of ^^Al but also 
affects the rate of its release in the ISM. Last, the coherent com- 
putation of the SNlb/c contribution (bv lLCOdh has led to an up- 
ward revision of the SN yields of the most massive stars. 

Nonetheless, when extrapolated to the subsolar metallicity 
that is indicated by some observations of the Cygnus complex, 
our predictions fail to account for the INTEGRAL/SPl measure- 
ments. Since no complete and coherent grid of nucleosynthesis 
yields for non-solar metallicity stellar models is available today, 
the discrepancy may well come from our extrapolation from the 
currently available data. One could also question, however, the 
metallicity estimate and the fact that a stellar population lying 
at about the same galactocentric distance as the Sun and being 



much younger could be subsolar. This discrepancy will not be 
solved until complete nucleosynthesis calculations at non-solar 
metallicities are performed and/or until the subsolar metallicity 
of the Cygnus complex is confirmed by a more recent estimate 
using improved stellar atmosphere models. These works are all 
the more expected because the solar metallicity has recently 
been revised from Z=0.02 to Z=0.014, which will obviously af- 
fect nucleosynthesis ( Asplund et al. 2005). In the meantime, our 
current understanding of massive stars and their nucleosynthe- 
sis will be tested further by applying the analysis presented here 
to other nearby OB associations with detectable gamma-ray line 
emission, like Sco-Cen. The potential effect of binarity on the 
predicted nucleosynthesis activity should also be evaluated. 

Finally, we performed a numerical simulation of the diffu- 
sion of ^*A1 inside a superbubble. Applied to the massive Cyg 
OB2 cluster, our model predicts that after 3 Myr, the ^^Al fills 
the entire superbubble and extends up to a radius of ~ 140 pc, 
in agreement with the observed extent of the 1809keV emission 
from Cygnus. Our simulation also indicates that most of the Al 
is affected by expansion velocities in the 50-300 km s"' range, 
which agrees with the velocities inferred from the 1809keV fine 
broadening. 
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